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Abstract. Discrete breathers are time-periodic, spatially localized solutions of the 
equations of motion for a system of classical degrees of freedom interacting on a 
lattice. Such solutions are investigated for a diatomic Fermi-Pasta-Ulam chain, i.e., 
a chain of alternate heavy and light masses coupled by anharmonic forces. For hard 
interaction potentials, discrete breathers in this model are known to exist either as 
"optic breathers" with frequencies above the optic band, or as "acoustic breathers" 
with frequencies in the gap between the acoustic and the optic band. In this paper, 
bifurcations between different types of discrete breathers are found numerically, with 
the mass ratio m and the breather frequency uj as bifurcation parameters. We identify 
a period tripling bifurcation around optic breathers, which leads to new breather 
solutions with frequencies in the gap, and a second local bifurcation around acoustic 
breathers. These results provide new breather solutions of the FPU system which 
interpolate between the classical acoustic and optic modes. The two bifurcation lines 
originate from a particular "corner" in parameter space (w, m). As parameters lie near 
this corner, we prove by means of a center manifold reduction that small amplitude 
solutions can be described by a four-dimensional reversible map. This allows us to 
derive formally a continuum limit differential equation which characterizes at leading 
order the numerically observed bifurcations. 



PACS numbers: 02.30.Oz, 05.45.-a, 45. 05. -fx, 63.20.Pw 



1. Introduction 

Discrete breathers are time-periodic, spatially localized solutions of the equations of 
motion for a system of classical degrees of freedom interacting on a lattice. They are 
also called intrinsically localized, in distinction, e.g., to Anderson localization triggered 
by disorder. A necessary condition for their existence is the nonlinearity of the equations 
of motion of the system, and the existence of discrete breathers has been proved 
rigorously for some classes of systems [Il[2],[3llll[5l[6l[7l[8]. In contrast to their 
analogues in continuous systems, the existence of breathers in discrete systems is a 
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generic phenomenon, which accounts for considerable interest in these objects from a 
physical point of view in the last decade. In fact, recent experiments could demonstrate 
the existence of discrete breathers in various real systems such as low-dimensional 
crystals [9], antiferromagnetic materials [10], Josephson junction arrays [llj, molecular 
crystals [I2], micromechanical cantilever arrays [13], and coupled optical waveguides 
[H]. The presence of discrete breathers has a significant influence on the properties of 
a system, and their relevance is discussed for a variety of physical effects, ranging from 
anomalous heat conduction to targeted energy transfer. For a review of some properties 
of discrete breathers see |15] . 

In the present paper, we study discrete breathers in the diatomic Fermi-Pasta- Ulam 
(FPU) chain, a chain of alternating light masses mi and heavy masses m2, coupled by 
anharmonic forces. This model has an acoustic and an optic phonon band, separated 
by a gap. For hard interaction potentials, discrete breathers are known to exist either 
as "optic breathers" with frequencies above the optic band, or as "acoustic breathers" 
with frequencies in the gap between the acoustic and the optic band [31 E] (existence 
results are also available for soft potentials [8]). In general, for some given parameter 
values (breather frequency, mass ratio, parameters of the interaction potential), a 
variety of different breather solutions can exist, and there have been attempts towards a 
classification of discrete breathers in certain models [16]. Here, we will study breathers of 
different types and bifurcations between them in the diatomic FPU chain, with breather 
frequency oo and mass ratio m = mi/m2 as bifurcation parameters. We will use both 
numerical and analytical methods to identify existence regions of different types of 
discrete breathers and their bifurcation lines in the small amplitude limit. 

In particular we numerically identify a period tripling bifurcation around optic 
breathers, which leads to new breather solutions with frequencies in the gap, and a 
second local bifurcation around acoustic breathers. These results provide new breather 
solutions of the FPU system which interpolate between the classical acoustic and optic 
modes. 

These numerical results are complemented by a local analytical study of the two 
bifurcations. One difficulty is that the two basic states consist in breather solutions, 
instead of having the FPU chain at rest as a starting state. Our approach uses 
the fact that the two bifurcation lines originate from a particular "corner" in the 
parameter space {uj'^^m). As parameters lie near this corner, we prove by means 
of a center manifold reduction that small amplitude solutions can be described by a 
four-dimensional reversible map. This allows us to derive formally a continuum limit 
differential equation which describes the numerically observed bifurcations at leading 
order. With this leading order analysis, the problem reduces to studying homoclinic 
bifurcations around reversible homoclinics of the asymptotic differential system. This 
allows one to compute the local shape of the "sector" in parameter space where 
breather bifurcations are numerically observed. In addition we provide leading order 
approximations of the new bifurcating breather solutions. 

After introducing the diatomic FPU chain in section [21 an overview regarding some 
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existence results for discrete breathers in this model is given in section [31 The first part 
of the paper is devoted to a numerical investigation, and we start by briefly reviewing the 
technique of numerical continuation of discrete breathers in section HI Numerical results 
obtained with this technique are presented in section [51 The second part of the paper 
contains the complementary analysis of small amplitude discrete breathers by means of 
analytic methods. To this purpose, the equations of motion of the diatomic FPU chain 
are reformulated as a map in some loop space in section 16.11 and a center manifold 
reduction, allowing to locally reduce the dimensionality of the map, is performed in 
section 16.21 Section [7] provides a formal asymptotic analysis for this reduced map. In 
section [HI this analytic result is confronted with the numerical findings, and we conclude 
with a summary and some comments in section [9l 

2. Diatomic FPU chain 

In 1955, E. Fermi, J. Pasta, and S. Ulam published a study on the failure of short- 
time energy equipartition observed for certain initial conditions in a one-dimensional 
Hamiltonian lattice with nonlinear nearest-neighbor interactions [I7j. Despite its simple 
form, this system reveals a surprisingly rich and complex behavior, and it has become 
one of the most studied models in the theory of nonlinear phenomena. Here we consider 
a diatomic extension of this model, consisting of a chain of alternate masses, coupled 
by identical nonlinear springs. The Hamiltonian function of the system is 



with momenta p„ G R and mass displacements g„ G R, and the interaction potential 
W is of the form 



where ooo, [3 G IR and O is the order symbol. We assume W to be an even function 
to simplify the computations, but all the analysis which follows could be extended to 
non even potentials containing a cubic term in the expansion (12|). Alternate masses 
correspond to the choices m2n+i = riii and m2n = ^2 for all n G where, without loss 
of generality, we choose mi < m2. Whenever resorting to numerical methods, we fix 
W (x) = + ^x^, mi = 1, m2 = (where m G (0, 1) is a parameter) and refer to 
a finite FPU chain with free boundary conditions, i.e., a Hamiltonian function 




,2 



(1) 




(2) 




(3) 



consisting of degrees of freedom. 



3. Discrete breathers in the diatomic FPU chain 



The first numerical observations of discrete breathers in diatomic FPU systems date 
back to the beginning of the 1990s [THl [TH]. A first rigorous proof of the existence 
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of discrete breathers in this system was obtained by Livi, Spicci, and MacKay [3j in 
1997. These authors showed for hard interaction potentials the existence of discrete 
breathers with frequencies above the optic band by continuation of a single-site breather, 
trivially existing in the "uncoupled" limit m = mi/m2 — >■ 0, towards small but non-zero 
mass ratios m. It was not until recently that the existence of small amplitude discrete 
breathers was proved for arbitrary mass ratios by James and Noble [3, [8] by means of 
a center manifold reduction. In the present article, this type of analysis is extended 
to the case of a higher dimensional center manifold. This extension allows to describe 
bifurcations between different types of discrete breathers occurring in certain parameter 
ranges. 

At least formally, a necessary condition for discrete breathers to exist appears to be 
the nonresonance of the discrete breather frequency uj with the phonon spectrum uOq of 
the linearized system, i.e., ujq/u} ^ N [20]. The plot of figure [H serves to illustrate where 
forbidden zones in parameter space (a;^,m), i.e., zones of frequencies uj not complying 
with this condition, lie. James and Noble proved the existence of discrete breathers 



m 




Figure 1. Bands (shaded) in parameter space (a;^,m) in which existence of discrete 
breathers is not expected due to resonance effects. The figure refers to the rescaled 
system ([6]), or equivalently to the original system ([T]) with mi = 1, rn2 — and 
= 1- The broad bands are the acoustic, respectively the optic, band of the linearized 
system. The narrower ones are images of the optic band, corresponding to breather 
frequencies lo whose multiples are resonant with the optic band. A countably infinite 
number of these images of the optic band exist. The images of the acoustic band are 
not plotted as they all lie inside the acoustic band and do not give any additional 
forbidden zones. Hatched areas sketch parameter values for which existence of discrete 
breathers has been proved for a hard interaction potential ('normally' hatched: proved 
in [3 [8], lightly hatched: proved in [3]). 
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outside (but very close to) the bands as introduced in the caption of figure [TJ For 
so-called hard potentials with (3 > 0, discrete breathers are found slightly above the 
bands, in the case of soft potentials with (3 < slightly below. For hard potentials the 
zones of proved breather existence are sketched in figured] (hatched areas). The result 
of Livi, Spied, and MacKay [3], in contrast, shows existence of discrete breathers for 
arbitrary uj"^ > 2, but close to the base line at m > (lightly hatched in figure [1]). 
Although rigorous existence proofs apply only to narrow regions in parameter space, 
numerical results indicate that discrete breathers can be found for arbitrary parameter 
values outside the forbidden zones [2T] . 

4. Numerical continuation of discrete breathers 

The focus of our interest lies on the various types of symmetric discrete breathers, their 
regions of existence in parameter space, and on the bifurcations occurring between them. 
In the first half of this article, with the intention to attain an overview of the types of 
discrete breathers that can be expected and on where bifurcations might show up, we 
resort to numerical techniques. The technique we use to compute discrete breathers 
up to numerical accuracy was put forward and discussed by Marin and Aubry in |22j . 
It is based on Newton's method to compute zeros of functions from a starting value 
sufficiently close to the desired zero. In this way, families of discrete breathers can be 
constructed by "continuing" a known breather solution of certain parameter value A [a 
vector, in our case (a;^,m), of any continuously varying parameters characterizing the 
solution] to a close-by value A'. This continuation procedure has been put on rigorous 
grounds by Sepulchre and MacKay [23j , by introducing the notion of "normal" periodic 
orbits which can be shown to persist for small continuous parameter changes. 

4.I. Seeds for numerical continuation 

With this continuation method at our disposal, and in order to compute a certain family 
of discrete breathers, this leaves us with the problem of finding a starting value (or 
"seed") for the numerical continuation, appropriate for the breather family of interest. 
The standard procedure to obtain such a seed is the anti-continuous limit, as mentioned 
in section [3] and utilized for the first existence proof [3] of discrete breathers in the 
diatomic FPU chain. In the limit m2 —>■ oo when every second mass is infinite, time 
periodic and spatially localized solutions exist trivially. The simplest one consists in the 
oscillation of one particle between "fixed walls" (two infinite masses), while all the other 
particles are at rest. However, for hard potentials ([2]) with (3 > 0, the continuation of 
such a seed will result in discrete breathers with frequencies above the optic band. This 
happens due to the fact that, in figured! such a seed corresponds to a point on the line 
(c(j^,m) = {uj'^,0) with u'^ > 2, and, as discrete breathers cannot be continued within 
the bands, a continuation from this seeds to frequencies with u"^ < 2 is in general not 
possible. 
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Figure 2. Snapshot of a gap breather at the instant of maxhnum oscillation 
amplitude for a diatomic FPU chain with coupling constants loq — I, [3 — 1 and 
(w^,m) « (0.256,0.126), computed numerically with a seed from the rotating wave 
approximation. The large masses m2 (large dots) oscillate with larger amplitude than 
the small masses mi (small dots). 



The bifurcations between different types of discrete breathers we are interested in 
will be found prevalently in the gap in parameter space (cj^, m) between the acoustic and 
the optic band. The seeds necessary to compute gap breathers by the above described 
numerical continuation method can be obtained by several strategies. 

One way is to use a rotating wave approximation which consists in approximating 
discrete breathers by using a small number of Fourier modes [2H [25l [26]. In the 
small amplitude limit, the slow spatial modulation of the fundamental harmonic is 
approximately determined by a nonlinear Schrodinger equation [271 [28]. This provides 
an explicit approximation of the envelope in the form of a soliton profile, which serves 
as an initial guess for the Newton method. For hard potentials, the computed starting 
states are of the "acoustic breather type", i.e., similarly to the acoustic modes of the 
linearized equations of motion, the large masses m2 oscillate at large amplitude, while 
the small ones mi move only slightly (see figure [2]). An entirely different strategy to 
obtain seeds for the continuation of gap breathers makes use of the anti-continuous limit 
with big masses m2 set to infinity in a slightly more sophisticated way than explained 
in section [H Instead of the simplest case with only one oscillating light mass mi, two or 
more oscillating ones with different frequencies are considered. Taking, for example, a 
light mass oscillating with period Ti, while its two adjacent light masses oscillate with 
T2 = |Ti, a discrete breather with period T = 3Ti is obtained (see figure [3] for an 
illustration). In this way, although Ti,T2 < ttV^, the overall period T of the discrete 

at rest — ^ ^ ^ » at rest 

Figure 3. Illustration of a seed for a gap breather from the anti-continuous limit with 
commensurate frequencies. The large balls indicate large masses m2 00, coupled 
to the small masses (smaller balls) by nonlinear springs. All masses except the three 
light ones indicated are at rest. The middle mass oscillates with larger amplitude, and 
therefore smaller period, than the two adjacent oscillating small masses. 
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breather can be larger than '^^^/2, giving rise to a gap breather. An example of a discrete 
breather produced by continuation of such a seed to finite mass ratios is provided in 
figure m In a similar way, more than three oscillating particles and other frequency 
ratios can be employed to obtain a large variety of gap breathers. 
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Figure 4. Snapshot of a gap breather at the instant of maxhxium oscillation 
amplitude for a diatomic FPU chain with coupling constants loq — I, (3 — 1 and 
(w^, m) « (0.253, 0.689), obtained by using a seed from the anti-continuous limit with 
commensurate frequencies 2uj and Sw. The large masses m2 (large dots) oscillate with 
smaller amplitude than the small masses mi (small dots). 

4-2. Boundary conditions and chain length 

As specified in section [21 free boundary conditions have been used for the numerical 
computation of discrete breathers. Under reasonable conditions on the chain length, 
i.e., on the number N of masses, this choice of boundary conditions does not affect 
the results reported: When the chain length is sufficiently large with respect to the 
localization length of the discrete breather, the oscillation amplitudes close to the 
boundaries are anyway smaller than the precision (typically of order 10~^^) of the 
numerical computation. As a consequence, the chain length necessary for the accurate 
numerical computation of a discrete breather increases with its localization length, and 
it is this need for large-A^-chains which imposes a limitation on the parameter values 
{uj"^, m) for which the numerical computation of small amplitude discrete breathers is 
possible. All numerical results plotted in this article were obtained for chains consisting 
of = 99 masses. 

5. Numerical results 

With the above described numerical methods at hand, bifurcations between various 
types of discrete breathers are observed and investigated: In section 15. ![ the 
bifurcation line of a bifurcation between symmetric and asymmetric discrete breathers 
is determined. In section 15.21 discrete breathers are investigated for parameter values 
(ci;^,m) close to a corner formed by the acoustic band and an image of the optic band 
(see figure [1]). Two kinds of bifurcations are detected in this region of parameter space: a 
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period-multiplying bifurcation on the one hand, and a local bifurcation around breathers 
whose oscillation is mainly performed by the large masses m2 on the other hand. Each 
of these bifurcations is of pitchfork type, and the bifurcation lines in the (tu^, m)-plane 
are computed in each case. All numerical results reported in this section are for coupling 
constants uq = 1, P = 1 in the interaction potential ([2]). 



5.1. Symmetric-asymmetric bifurcation 

At the end of the 1990s, a bifurcation between symmetric and asymmetric discrete 
breathers in the diatomic FPU chain was reported by Cretegny, Livi, and Spicci |21j . 
These authors investigated two relatively simple types of discrete breathers, the so- 
called Sievers-Takeno mode and the Page mode (see figure [5] for an illustration of these 
modes), above the optic band. They performed a linear stability analysis of these modes. 




Figure 5. Snapshots of breathers from the same breather family at the instant of 
maximum oscillation amplitude for a diatomic FPU chain with coupling constants 
cjo = 1, /3 = 1 for parameter values uj^ = 5.8397 and m e {1.0000,0.9920,0.9886}. 
The linearly stable antisymmetric, bond-centered breather (Page mode, left figure) at 
mass ratio m = 1 deforms into an asymmetric breather for m < 1 (middle figure) . For 
sufficiently small mass ratios, this branch of asymmetric breathers and the branch of 
symmetric, site-centered breathers (Sievers-Takeno mode, right figure) collide and a 
bifurcation takes place. 



observing that for small values of the mass ratio m the Sievers-Takeno mode is linearly 
stable, while for m = 1 the Sievers-Takeno mode is unstable and the Page mode is 
linearly stable (see [21] for details of the stability analysis). This observation led to the 
conjecture, and the subsequent numerical confirmation, of the existence of a bifurcation 
connecting these two branches of solutions. We give more informations about this 
bifurcation by computing numerically the bifurcation line in the (cu^, m)-plane. 

Since the bifurcation takes place above the optic band, the anticontinuous limit 
can be used as a seed for the numerical continuation. Discrete breathers of the Sievers- 
Takeno-type are obtained from the "standard" anticontinuous limit of alternate finite 
masses mi and infinite masses m2, where only one of the light masses mi oscillates 
while all the others are at rest. Numerically, the continuation is possible across the 
entire range of mass ratios ^ m ^ 1. 

The computation of the Page mode in the diatomic chain is slightly more complex 
using the anticontinuous limit. This requires a preliminary step, where discrete breathers 
of the Page-type are produced from an anticontinuous limit where every third mass is 
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Figure 6. Bifurcation line in the parameter plane (uj'^,rrij for the symmetric- 
asymmetric bifurcation of discrete breathers. Below the bifurcation line, linearly 
stable symmetric site-centered breathers (Sievers-Takeno mode) are observed. At 
the bifurcation line these breathers become unstable and, additionally, linearly stable 
asymmetric breathers show up. For equal masses m = 1 these asymmetric ones 
transform into anti-symmetric bond-centered breathers (Page mode) . Lines connecting 
the numerically determined data points are drawn to guide the eye. 

infinite while all other masses are equal and finite. Continuing a seed, with two adjacent 
light masses oscillating out of phase and all others at rest, up to mass ratio m = 1 (i. e., 
all masses in the chain are equal), the Page mode is obtained. Then, starting from mass 
ratio one, the Page mode can be continued to the diatomic chain of alternate masses. 

In figure [6l the bifurcation line for the bifurcation between symmetric and 
asymmetric discrete breathers is plotted in a section of the (cj^, m)-plane. Below the 
bifurcation line, discrete breathers of the Sievers-Takeno type are observed and are found 
to be linearly stable. At the bifurcation line, these breathers loose linear stability and, as 
additional solutions, asymmetric breathers emerge. Above the bifurcation line, unstable 
Sievers-Takeno type breathers and linearly stable asymmetric breathers coexist. At mass 
ratio m = 1, the asymmetric breathers transform into antisymmetric breathers of the 
Page type. Note that, in the diagram sketching the phonon bands of the linearized 
equations of motion (figure [1]), the bifurcation line emanates from the endpoint of the 
upper band edge of the optic band. It is a common feature of all the bifurcation lines 
investigated in this article that they emanate from endpoints of bands or corners formed 
by bands and their images as shown in figure [U 

5.2. Bifurcations close to "corners" formed by bands and their images 

In reference [TJ IH], James and Noble proved the existence of discrete breathers for 
the diatomic FPU chain in the vicinity of edges of the phonon bands, with different 
functional forms of the discrete breathers close to the different band edges. At the 
"corners" in parameter space (cj^, m) formed by the acoustic band and the images of 
the optic band (see figured]), one might conjecture (for hard potentials) the occurrence 
of bifurcations connecting these two types of breathers, i.e., an optic and an acoustic 
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breather. This conjecture was the original motivation for the study of gap breathers 
reported in this article. 

We chose the corner formed by the acoustic band and the "second image" of the 
optic band at (a;^,m) = (|, |) (see figure [1]) for an exemplary investigation. For 
the numerical detection of bifurcations in this region we utilized seeds obtained from 
the rotating wave approximation for the numerical continuation. Two types of local 
bifurcations of pitchfork type were found: 

(i) a period multiplying bifurcation (local bifurcation around an optic breather whose 
oscillation is mainly performed by the small masses mi), 

(ii) a local bifurcation around an acoustic breather whose oscillation is mainly 
performed by the large masses m2- 

As parameters are varied, these two bifurcations connect discrete breathers with 
prevalent oscillations of the heavy masses m2 to discrete breathers with prevalent 
oscillations of the light masses mi. In the left part of figure [71 the transformations of the 
shape of discrete breathers successively undergoing the two bifurcations are illustrated 
for two different bifurcation branches. The right column of figure [7| compares these 
numerical solutions to analytical approximations obtained in sections [7] and [HI (these 
results will be commented in section [8]). 

As a measure of how strongly the oscillation of the discrete breather is dominated 
by the heavy or the light masses, the quantity 

^ sup J2 lL (t) 

^ = o In * " 2 TIT (4) 

2 sup E g2„+i (t) 

t n 

may be used. We consider here the ratio r of the mean-square displacement of heavy 
masses by the mean-square displacement of the light ones, whose logarithm is given 
hj g = Inr. A positive value of g corresponds to a discrete breather dominated by the 
oscillation of the large masses, whereas a negative value of g indicates that the oscillation 
of the small masses is dominating. Plotting g for a family of discrete breathers (i.e., 
as a function of the parameter labeling the members of the family), the interpolation 
between acoustic and optic breathers provided by the bifurcating solutions can be nicely 
illustrated (see figure [8]). We have plotted g instead of r — 1 in order to stretch the region 
where r < 1, thereby increasing the visibility of this region in the plot. 

For both, the local bifurcation around an acoustic breather and the period tripling 
bifurcation, the bifurcation lines in the parameter plane (u;^,m) have been computed 
numerically and are plotted in figure M Rehable data points even closer to the bands 
could not be obtained numerically, as the localization length of the discrete breather is 
increasing when approaching the band and therefore size effects in the finite systems 
used for the numerics become non-negligible. 
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Figure 7. Snapshots of gap breathers from the same breather family at the instant 
of maximum oscillation amplitude for a diatomic FPU chain with coupling constants 
tjQ = I7 /3 = 1 for different points on the line (w^, m) = (2m + 0.01555, to) in parameter 
space. For to « 0.1428, the large masses TO2 (large dots) have a larger oscillation 
amplitude, while for decreasing m the amplitude of the small masses' oscillations (small 
dots) increases. The left and middle columns are the results of a numerical computation 
of discrete breathers from two different bifurcation branches. The right column shows 
approximate analytic results as obtained for the indicated parameter values from 
equations ([50|) . ((55|l . ([5T|l . and (|48|) (from top to bottom) for comparison with the 
numerical results in the middle column (see section [8] for details) . The qualitative 
agreement of numerical and analytic results is excellent throughout. The quantitative 
agreement, as expected, is good only for weakly localized breathers (e. g. in the upper 
row) for which the assumptions of the analytic calculation are met. For parameter 
values (w^, to) even closer to the phonon band edges for which also the breather types 
shown in the lower row would be less localized, we encountered problems with the 
convergence of the algorithm used for the numerical computation and therefore could 
not compare analytic and numerical results. 



6. Center manifold reduction near critical parameter values 

The bifurcations observed numerically in the previous section were found to lie in the 
vicinity of crossings of phonon band edges, supporting the conjecture that bifurcation 
lines emanate from the corners formed by these bands. In this section we apply the 
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Figure 8. Along the line (u;^,m) = (2m + 0.01555, to) in parameter space, the 
quantity g defined in (jj]) is plotted for discrete breathers from different bifurcation 
branches as a measure of how big the contributions of the light and the heavy masses, 
respectively, are to the overall oscillation. The numerical continuation of solutions 
connects discrete breathers with predominant heavy and light mass oscillations (top 
and bottom lines) via a succession of two pitchfork bifurcations. Computations are 
performed for a diatomic FPU chain with = 99 degrees of freedom and coupling 
constants loq — 1 and /3 — I. 



technique of center manifold reduction for an analytic investigation of these bifurcations. 
When proving the existence of discrete breathers close to the phonon band edges in the 
diatomic FPU chain in [71 [8], James and Noble rewrote the equations of motion of the 




Figure 9. Bifurcation lines close to the corner at (a-'^,TO) — (|, |) for the period 
tripling bifurcation (lower lines and 0) and the local bifurcation around an acoustic 
breather (upper lines and -I-). Above the upper bifurcation line, discrete breathers with 
predominant heavy mass oscillations are observed. At the bifurcation line, additional 
breathers show up with light mass oscillations growing in importance. At the lower 
bifurcation line, the frequencies of the breathers with light mass oscillation are tripled 
which, in the plotted parameter plane, corresponds to a jump to frequencies above the 
optic band (see figure [1]). Dashed lines connecting the numerically determined data 
points are drawn to guide the eye. The solid lines arc the analytic approximations 
pO)) (lower line) and (|47)) (upper line) of the bifurcation lines obtained by means of a 
center manifold reduction, valid in the vicinity of the corner. 
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diatomic FPU chain as a map in some loop space and accomplished a reduction of the 
spatial dynamics onto a two-dimensional center manifold. In the vicinity of the "corners" 
in parameter space we are interested in, a similar analysis can be performed, with the 
complication that the relevant center manifold will be four dimensional. In section [7] we 
shall analyze asymptotically the four-dimensional map on the center manifold. A set 
of equations will be obtained whose solutions correspond to different types of discrete 
breathers, some of which can be identified as the ones observed numerically in the 
previous section. 

6.1. Formulation of the FPU system as a map in a loop space 

We consider the diatomic FPU chain ([T]) consisting of an infinite number of degrees 
of freedom, and reformulate this system as a map in loop space. We start from the 
equations of motion 

d2 



rrir. 



dt^ 



qn = W (g„+i - g„) - W (g„ - g„_i) , n G Z, (5) 



where m2n+i = = ^2 {f^i < ""^2)? and the g„ represent mass displacements 

from their reference position. W is the interaction potential as defined in ([2]), and a 
prime denotes its derivative. In order to reduce the number of parameters, we rescale 
([5]) by setting g„ (t) = UqX^ [uot/^ym^), which yields 



(6) 



-1 / 

TTT' -^^2n = V (-^2n+l " X2n) — V (-^2n " ^2n-l) , 
d2 

-^X2n+1 = V {X2n+2 — -^2n+l) — V {X2n+1 — -^2n) , 



with mass ratio m = mi/m2 G (0, 1) and 

V{x) = ^x^ + f^x^ + 0{x^). (7) 

Being interested in discrete breathers, we look for time-periodic solutions of ([6]) with 
frequency uj. In order to work in a space of functions with fixed period 2it, we rescale 
Xn by setting Xn (t) = Xn (cut) (x„ being 27r-periodic), which transforms the equations 
of motion ([6]) into 

d^ 
-1, .2 ^ 



m LU 



dt 



X2n — y' ix2n+l — a^2n) ~ V' ix2n — X2n-l) 



(8) 



;X2n+l — V ix2n+2 ~ X2n+l) — V {X2n+1 — X2r. 



dt^' 

Breather solutions of ([6]) with frequency uj correspond to 2'7r-periodic solutions of ([8]) 
satisfying the localization condition lim|„|^oo x„ (t) = 0. For the center manifold 
reduction procedure, however, we do not need assumptions on the spatial behaviour 
of solutions, except the assumption of small amphtude. 

Using the fact that V is locally invertible, one can formulate dHj) as a first order 
recurrence relation in a loop space. Indeed, the first equation in ([8]) defines a map 
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(X2n, X2n-l) 

consider u„ 



with 



X2n+i, the second one a map 

3^2n; "^n = X2n-1 ^nd = {Un, Vn 
Fm,u} {y^n) 



[x2n,X2n+l) ^ X2n+2- MoiC CXphcitly, 

. Then, (IHl) can be rewritten as 



and 



G^2 [G^2/^ {u,v) ,u] 



v 



u 



v) + a 



d^ 

dt2 



+ U. 



(9) 



(10) 



:iii 



Note that the fixed point 



of Fm,Lo corresponds to the FPU chain at rest. In 



the following, m and uo play the role of bifurcation parameters. Even if the recurrence 
relation ([9]) is infinite-dimensional and ill-posed (the operator -F^,^ is defined on a 
domain D specified below), small amplitude solutions lie on a finite-dimensional center 
manifold whose dimension depends on the number of resonant phonons in the parameter 
regime considered [HI [H [7]. The reduction to this finite-dimensional center manifold is 
the crucial step which allows for an analytic treatment of discrete breathers close to 
certain phonon band edges. 

As a next step, we define appropriate function spaces on which Frn,u] is acting. We 
restrict our attention to even functions of t, thereby reducing the dimension of the center 
manifold by a factor of two (half of the Fourier modes are eliminated). Moreover, we look 
for 27r-periodic oscillations consisting entirely of odd order harmonics cos [{2k + 1) t] with 
/c G No, which is possible due to the evenness of V . To be more precise, we introduce 
the function spaces (n ^ 0), defined by 



HI 



{u e /f" {R/27tZ) I u (-t) =u(t) ,u(t + n) = -u (t)} 



(12) 



where if" denotes the standard Sobolev space. For Fm^^ {u, v) to be well-defined, we 
look for (u„, Vn) G D = x H^, and then the recurrence relation holds in 
X = X H^. Due to the evenness of V the operator Fm,uj maps D onto X (even 
order Fourier harmonics are not generated by applying F^^^), and Fm,uj is analytic in a 
neighborhood of {u,v) = in D. 

For even potentials, working in a function space with only odd order Fourier 
harmonics allows to extend the existence results obtained by James and Noble for 
acoustic breather solutions (see [8], theorem 1 iii)). This procedure suppresses half 
of the forbidden zones in figure [Tj namely the ones corresponding to even multiples 
of the optic band. In the general case of asymmetric potentials, where all Fourier 
harmonics are involved, the map F^^u) linearized at (m, v) = possesses a pair of 
imaginary eigenvalues on the unit circle when parameters lie in these forbidden zones, 
which causes an obstruction for a breather existence result (see [8]). These additional 
eigenvalues (corresponding to even order Fourier harmonics) are eliminated with the 
present functional setting. 
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6.2. Reduction of the dimension near a resonance of order three 

As for the numerical investigation in section 15.21 we choose for a detailed analysis as 
an example the region of the crossing of the upper band edge of the acoustic band 



frequency resonance of order three between the optic phonon with wave number q = 
and the acoustic phonon with wave number q = tt. 

In order to study the local dynamics of ([9]) around F„ = 0, we first need to determine 
the spectrum of the linearized map DF^^w (0) close to the unit circle. This linear analysis 
has been already performed in reference [H] using a different (but completely equivalent) 
reformulation of the FPU system. In what follows, we briefly recall the useful spectral 
properties in our analysis and refer to reference [8] for more details. For all parameter 
values (m, u), the spectrum of DF^^^ (0) on the unit circle consists in a finite number of 
eigenvalues which can be deduced from the usual dispersion relation for linear phonons 
(see below). Moreover, this central part of the spectrum is isolated from the hyperbolic 
part. Therefore DF^^uj (0) satisfies the property of spectral separation required to apply 
the center manifold theorem [6]. 

D-Pm,aj (0) admits a pair of simple eigenvalues e^*'' on the unit circle (with an 
associated invariant subspace spanned by the vectors (cos (kt) , 0) and (0, cos (kt)), for an 
odd integer k) if ku lies inside one of the acoustic or optic phonon bands plotted in figure 
[H With this pair of eigenvalues, equation ([6]) linearized at X„ = admits solutions of 
the form {X2n, -^2n+i) = cos (kut) [e*^" ( + e~*^" C] for appropriate eigenvectors ( G C^. 
These solutions correspond to linear standing waves with wave number q, where q and 
ku) need to satisfy the classical dispersion relation for the diatomic FPU chain. This 
observation allows us to determine the structure of the spectrum in a rather simple 
way. For {u!,m) = (|, |), the part of the spectrum of DFm^ui (0) which lies on the unit 
circle consists only of two double non-semi-simple eigenvalues —1 and +1, which split 
generically into two pairs of simple eigenvalues as parameters are slightly varied. The 
invariant subspace associated with eigenvalue —1 is spanned by the vectors (cos (t) , 0) 
and (O,cos(t)), whereas the invariant subspace belonging to eigenvalue +1 is spanned 
by (cos (3t) , 0) and (0,cos(3t)). If {uj,m) = (|,|), uj lies at the top of the acoustic 
phonon band (corresponding to the wave number g = vr), Scu lies at the top of the optic 
band (wave number g = 0), and no higher multiple of u lie in the phonon spectrum. 
With this choice of parameters (a;,m), the central subspace (i.e., the subspace which 
is invariant under DFm,uj (0) and which is associated with the eigenvalues lying on the 
unit circle) is four-dimensional. 

In the following, we set 



with /i, 77 ^ 0, which implies u; ~ | and m ^ |. Equation ([9]) [see also ([H])] can be 




m ^uj'^ = 2 + /i. 



(13) 
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rewritten as 

d2 



d 



(14) 



dt 

The center manifold theorem |6] which has been apphed in [TIE] states that, for yU, r/ ^ 0, 
small amplitude solutions of f[T^ (with (m„, i;„) G D for all G Z) have the form 

Un = an COS (t) + On COS (3t) + [(f (a„, c„, bn, /i, r])] (t) , 

(15) 

f„ = dn COS (t) + 6„ COS (3t) + (a„, d^, c„, 6„, /x, r^)] (t) . 

The functions (v?, ip) : ^ D are combinations of higher order Fourier components 
cos((2A; + l)t) with k ^ 2, they depend smoothly on = {an,dn,Cn,bn), fJ^y^i, and 
contain only higher order terms in Y^,fi,7], i.e., 

V {0,0,0,0, fi,r]) = 0, Dyc(^(0,0,0,0,/i,r/) = (16) 

(the same holds for ip). A four-dimensional center manifold, locally invariant under Fm,u) 
and containing all small amplitude solutions, is locally constructed as the graph of the 
function {ip,4') [■,^,ri] defined on the central subspace. As indicated by (fT6|l . the center 
manifold is tangent to the central subspace at F„ = (even for (/i, ?]) ^ 0, since the 
harmonics cos (t) and cos (3t) span an invariant subspace of the linearized dynamics for 
all parameter values). 

The center manifold reduction theorem can be intuitively understood by considering 
the linearized system at F„ = 0, given by F„+i = DF^^^ (0) Yn for (a;,m) = (i, |). In 
the linear case, all bounded solutions lie on the 4-dimensional invariant central subspace 
(all the other solutions diverge exponentially as n +oo or —oo since they possess 
a component along an hyperbolic mode). A similar result holds true locally in the 
nonlinear case, where the central subspace is replaced by a 4-dimensional invariant local 
center manifold containing all small amplitude solutions. 

We note that ip (■, fi, rj) and ip fi, rj) commute with — / (the symbol / denotes the 
identity operator) since the center manifold is invariant under the symmetry — / of ([9]), 
which originates from the evenness of V. Consequently, ip and ip are of order Hl^'^H^ as 
— > in R^, where || ■ || denotes the Euclidean norm. In the following, we denote by 
= {(f, Ip) the reduction function having the center manifold as its graph. Furthermore, 
we denote by Pc the spectral projection on the central subspace, 

+ 00 

Pc{^ Ak COS {{2k + I) t)^ = Ao cos (t) + Ai cos (3t) , (17) 

fc=0 

and denote by X the isomorphism 

X {a, d, a, bY = {a cos {t) + ccos (3t) , dcos (t) + 6 cos (3^))""" (18) 

between M.'^ and the central subspace. For ^,ri ^ 0, the sequences of coordinates Y^ G 
of small amplitude solutions on the center manifold are given by the four- dimensional 
recurrence relation 

^ra+l — ffJ.,V O^n) 5 (19) 
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where 

Z^,, (n = j-ip, F„,.(xy'= + (F^ fi, v)) (20) 

for all G M"^ in a neighborhood of 0, and where the parameters m, uj are defined 
by ( |T3i) . Equation (fT9l) is obtained by substituting Ansatz ( |T5l) into (IHl) and projecting 
the resulting equation onto the central subspace. Note that (fT9!) inherits the invariance 
of ([9]) under — /. Moreover, the invariance of ( JT^ under the symmetry (u„, f„) — »• 
(-M_„, -f-„+i) induces the invariance (a„, c„, 6„) (-a_„, -c/_„+i, -c_„, -6_„+i) 
in the recurrence relation ( fTOl) . Since = O (||y^||^) as ^ 0, we have in addition 

z^,, (n = j-ip, (xn + o (iiri^) . (21) 

In what follows we compute the principal part of the reduced mapping ( |T9l) . Due to 
the complicated structure of -Fm.oj, we shall not directly compute the Taylor expansion 
of (12T1) and proceed instead in the following way. Substituting Ansatz (fT5l) into (fl^ . 
projecting onto the central subspace (i. e., applying P^) and identifying the components 
along the harmonics cos (t) and cos (3t) results in the following recurrence relations 
(computations are lengthy but straightforward) 

dn+l +dn + flttn = -Jp {dn - an) [(rf„ - a„)^ + ((i„ - a„) {bn - C„) + 2 {bn - Cnf] 

-|/3 {dn+l - an) [{dn+l " anf + ((in+1 " On) (6^+1 - C„) + 2 (6„+i - Cnf] 

+0((||V:il + ||i;VII)'), (22a) 

&n+l + &n + (9/U + 16) C„ = -(3 [\ {dn - anf + | (rf„ - anf {bn - C„) + | (6„ - Cnf] 
-(3 [\ {dn+l - anf + I {dn+l " anf {bn+1 - C„) + | (6„+i - Cnf] 

+o((i|i;ii + i|i;vii)'), (226) 

an+1 + an + (77 - I) (i„+i = |/5 ((i„+i - a„+i) 

X [(c^n+l — an+lf + ((in+l " C^n+l) {bn+1 — Cn+l) + 2 {bn+1 — Cn+lf] 

+ (c^n+1 — O^n) \{dn+l — anf + {dn+l ~ Q-n) (&n+l " C^,) + 2 (6„+i — C„)^] 

+0((||K„1| + ||i;VII)')' (22c) 
Cn+i + c„ + (i + 9?7) 6„+i = 

/5 {dn+l - anf + I {dn+l " a„)^ (6„+i - c„) + | (6„+i - Cnf] 

+13 ((in+l ~ O-n+lf + f (c^n+1 ~ O-n+lf {^n+l C„+i) + | (6n+l ~ C„+i)^] 

+0((||K„1| + ||i;VII)')- M 
Equations fl22ap - fl22dj) represent an implicit form of the recurrence relation f[T^ . Using 
the implicit function theorem, one can reformulate (I22ap - (]22dj) locally as a map 

//^,r; • {O'n, dn, Cn,bn) ^ {an+1, dn+l, Cn+l, bn+l) ■ 

In conclusion, in the parameter regime /i, r/ ~ 0, the center manifold theorem 
from [6] guarantees that small amplitude solutions of (HM are determined by a 
four- dimensional mapping whose principal part is given by (l22^ -( l22d|l (in addition 
the reduced mapping preserves the symmetries of the infinite-dimensional system). 
Consequently one can state the following theorem. 
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Theorem 1 There exist neighbourhoods Q, U and A of {u, v) = in 3, Y"^ = in 
and {fi, rj) = in R^, respectively, and a C^-map (p : ^ D with (p (0, 0, 0, 0, /i, rf) = 0, 
Dyc0 (0, 0, 0, 0, fi,r]) = such that the following holds for all {fi, rj) G A; 

(i) All solutions of [I4\ ) such that {un, fn) G ^ for alln eJj have the form ( 173]) . where 
= {an, dn, Cn, bn) satisfics a 4- dimensional recurrence relation with principal part 
given by ^MS^-WM)- 

(a) IfY^ = (a„, dn, Cn, bn) is a solution of Ii22a\) - (22d\) with Y^ for all n E Z, then 

l[T5\) defines a solution of [14\ )- 
(Hi) The invariance of under the symmetries {un, Vn) — > {—Un, —Vn) and {un, f„) — >■ 

{-u-n, -v-n+i) induces the invariance under Y" — > -Y" and {an, d 

n , Cji , On ) ^ 

(— a_„, —d-n+i, —C-n, — &-„+i) for the reduced system Ii22a\) - (22d\) . 

Small amplitude orbits of the reduced map homoclinic to correspond via theorem 
[T]to exact discrete breather solutions of the FPU system. In what follows we shall not 
prove the existence of such exact homoclinic orbits. Instead we shall derive in the next 
section a continuum limit differential system which provides approximate homoclinic 
orbits of the map. These analytical results will be used in section [8] to explain the 
numerically observed breather bifurcations. 

In addition, the leading order continuum limit of the reduced map obtained in the 
next section might be used, together with reversibility [theorem [1] (iii)], as a tool for 
proving the existence of exact small amplitude homoclinics close to the continuum limit. 



7. Continuum limit and homoclinic bifurcations for the reduced map 
7.1. Continuum limit of the reduced map 

We assume, as before, that fj,,ri ^ with /x > 0, and introduce a small parameter e, 
setting /i = and rj = ne^. For e = 0, the recurrence relation fl22al) - (l22dl) linearized at 



the fixed point (a, rf, c, 6) = reads 

= dn+i + dn, (23a) 

= + 6„ + 16c„, (236) 

= an+l + Cln - i^n+l, (23 c) 

Q = Cn+l + Cn + \bn+l, (23(i) 

with bounded solutions given by 

{an, dn, Cn, bnV = « (("I)" , 0, 0, 0)^ + b (O, 0, -|, l)"" , (24) 



a, 6 G M being arbitrary. Indeed, for e = the spectrum of DFm.uj (0) on the unit circle 
consists in a double non-semi-simple eigenvalue —1 [with the corresponding invariant 
subspace spanned by cos {t)] and a double non-semi-simple eigenvalue [invariant 
subspace spanned by cos (3t)]. 
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We search for small amplitude approximate solutions of fl22ap - (122dj) in the form of 
a slow modulation of fl24l) supplemented by higher harmonics, i.e., 

/ a„ \ /(-ir\ / \ 



v 













((-!)" + (25) 



k>2 



V 1 / 



dn+1 

where ^ = en. This approach has some analogy with the multiple scales expansion 
technique (see j2S] and its references), with the fundamental difference that the FPU 
system has been locally reduced for time-periodic solutions (by theorem [1]) to a finite- 
dimensional problem. In what follows we compute the principal terms in the expansion 
fl25|) . Inserting fl25|) into the recurrence relations fl22ap - fl22d|) and identifying the different 
harmonics and powers of e, we obtain 

a„= (-l)"eAo(en) + 0(e2) 
e5o (en) + 0(e2). 



dn 



leBo (en) + 0(e2), 
-l)"e^^(6n) + 0(e^) 



where Aq, Bq are solutions of system 
d'Ao 











de 



Ao + f(3 {A', + §AoB',) 



- 9 (i + 16^) B, + {2AlBo + ii?o') . 



(26a) 
(265) 
(26c) 

(26d) 

(27a) 
(275) 



Details of this calculation are reported in the appendix. In this way we have obtained 
a formal continuum limit f l27ap -( l275j) . valid for fi,r] ^ 0, of the four- dimensional 
recurrence relation fl22al) - (]22dj) . The corresponding approximate solutions (I26al) - (]26d[l 
of the reduced recurrence relation in turn correspond to approximate solutions of the 
FPU system f|T^ having the form 



;-!)" eAo {en) cos (t) - ^e^o (en) cos (3t) + 0(e^ 



(28) 



eBo (en) cos (3t) + O(e^). 

Normalizing the solutions of the system f l27ap -( l275l) by setting -Bo (0 = (0 ^^d 
^0 (0 = (0 5 the equations read 
d^a 

■^-la + ^ha^ + Qhab'^ = 0, (29a) 

— -9(^ + 16K)b + 6ha% + ^hb^ = 0, (29b) 
at, \4 / I 



with 



(30) 
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In this notation, approximate solutions, valid for /i, ^ 0, of the original equations 
of motion of the diatomic FPU chain are given by 

(t) = (—1)" (e^) COS (tot) j=eb (m) cos (Scot) + O(e^), 

(31) 

X2n-i (t) = ej^b (en) cos {3ut) + 0{e'). 

The small parameter e (which measures the amplitude of the solutions), the frequency 
00, the mass ratio m and the parameter k in (12961) are linked by the relations 

^2 

— = 2 + = T + «^e^- (32) 

m 4 

The asymptotic analysis leading to this result is a priori valid for e ~ 0, k being fixed. 

This corresponds to the regime u ^ ^, m ^ ^, with m/cu^ < |. 

7.^. Homoclinic bifurcations in the continuum limit 

In this section we analyze some local homoclinic bifurcations for system fl29al) - (]296p . 
These bifurcations show up for hard interaction potentials, i.e., we shall assume /3 > 
in the potential in ([7]). 

Homoclinic bifurcations in systems of the type (I29al) - fl296p have been analyzed by 
several authors in different contexts (see e.g. [301 Ell [32] ) . For example, such a system 
can be obtained when considering two coupled nonlinear Schrodinger equations (CNLS 
system) and one is looking for time-periodic solution components. The CNLS system 
has been introduced to describe coupled mode dynamics in different physical contexts, 
such as fluid dynamics or nonlinear optics (see e.g. references in [3lj). In the context of 
nonlinear lattices, the CNLS system has been derived in reference ^7] for the Kac-Baker 
model, in order to describe the nonlinear modulation of two superposed linear modes 
with equal group velocities (this analysis is based on a formal multiscale expansion). 

Here we consider pitchfork bifurcations from a single-component homoclinic 
orbit of f l29al) -( l29B ((a, 0) or (0,6)) to two-components homoclinic orbits (a, 6). 
These bifurcations have been analyzed by Yang [30] using a classical perturbative 
analysis for pitchfork bifurcations (see e.g. [33]), and the computations in reference 
[30] can be readily employed here. As will be shown in the next section, these 
bifurcations correspond (qualitatively and quantitatively) to the numerically found 
breather bifurcations in section 15.21 



7.2.1. Homoclinic orbits bifurcating from {0,b): We first normalize the coefficients of 
f l29ap -( l296i) . We assume k > —-^ and set 

6(0 = aoM(50, a{0 = (^ov{SO, (33) 

where Oq = [ 27/1 (i ^^f^) ]^^^ 6 = 3 + IGk)^^^. This yields 

d^u -? ^ 9 

-u + w^ + /Squv^ = 0, 

(34) 



dx^ 



da;2 



- Q^v + pou^v + = 0, 
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where x = 6^, f3o = ^, 7 = (^)^ and = ^ (i + IQk) System has the 
homochnic solution (m, v) = {uh, 0) with Uh (x) = sech (x). We shall look for critical 
values of fl at which solutions homoclinic to bifurcate from {u,v) = {uh,0). 

This problem has been previously analyzed by Yang [30j for a slightly more 
symmetric system (with 7 = 1). This analysis readily applies here (with only a few 
straightforward modifications of the coefficients) and we refer the reader to [3U|| for 
details (see also [32] for more details on the linearized problem at {u,v) = {uh,0)). 
There exists a finite number of critical values 

Qn = s — n, (35) 

where s = (a/1 + SPo — l) and n ^ is an integer. A pitchfork bifurcation occurs 
from {u,v) = (m/i,0) at any of these critical values, due to the invariance v — > —v of 
f l29ap -( l296|) . If n is even, bifurcating solutions are even in x. If n is odd, u is even 



and V is odd for bifurcating solutions. The number of critical values depends on the 
size of Po, and is equal to one in our case. We have Qq ~ 0.272 which corresponds to 
K = Ko^ 0.149. 

The branch of homoclinics bifurcating from {u, v) = {uh, 0) as $7 ~ Qq can be locally 
parameterized by 

u = y2sech (x) + O (A^) , (36) 
V = Asech^(x) + (|A|=^) , (37) 

= nl + 9x^ + {x''), (38) 

where A is a small parameter (see [30j, paragraph 3.1). The sign of 

sech^"* (x) dx 

determines whether the pitchfork bifurcation is sub- or supercritical. Here 7— ~ —0.01 
and consequently one has < fig and k, > Kq for bifurcating homoclinic solutions. 
Returning to the original parameters {uj'^,m) given by fl32l) and fixing u"^ > 1/4, 
bifurcating homoclinics exist for 

m > mpT (^) = I + «PT (w' - i) + O [{lu^ - i)'] , (40) 

where = ^ — ~ 0.08, and the index PT is for "period tripling". 

7.2.2. Homoclinic orbits bifurcating from {a,0): We again normalize the coefficients of 
f l29ap -( l296j) . We assume k > — ^ and set 



a(0 = aiw(^e), &(0 = «i^(^e), (41) 
where Oi = ( |^ )^^^- This yields 



d^u 
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where a; = = , 7 = (f )^ and = ^ (i + System (02]) is of the same 

type as ( IMl) . hence the analysis performed in the previous section still applies. We look 
for critical values of Q at which solutions homoclinic to bifurcate from {u, v) = {uh, 0). 
These critical values are given by 

Qn = 's — n, (43) 



where s = | ( \/ 1 + 8/3o — 1 ) ~ 6.322 and n is an integer with ^ n ^ 6. 



In what follows, we concentrate on the bifurcation near Qq = 's, which corresponds 
to K = ~ 0.47. The branch of homoclinics bifurcating from {u, v) = {uh, 0) as ~ Qq 
can be locally parameterized by 

u = y2sech(a;) + (A^) , (44) 
V = Asech^(a;) + (|A|^) , (45) 

n^ = hl + ex^ + o (a^) , (46) 

where A is a small parameter. Since 6 has the same sign as 7 — s"^, one has < 
and < i^o for bifurcating homoclinic solutions. Therefore k, < on the bifurcating 
solution branch. Returning to the original parameters (c(j^,m) given by (132|) and fixing 
tu^ > 1/4, bifurcating homoclinics exist for 

m<muMiuj) = l + auM{uJ^-\)+0 {io^ - \f , (47) 

where auu = \ — ~ 0.367, and the index HM is for "heavy mass oscillations". 

We end with a remark concerning the above pitchfork bifurcations in fl29al) - (]296p . 
related with the invariances (a, h) (—a, h) and (a, h) (a, —h). This system formally 
corresponds to ([22aj)-([22dj) in a continuum limit, but the corresponding symmetries 

iflni dn, Cn, bn) ^ ( dny d^, On, bn) and (ctji, d^, On, bn) ^ (^n, Cn, &n) are lost 

for the exact system fl22ap - fl22dj) . Consequently, we only expect imperfect pitchfork 
bifurcations (i.e., with nonsymmetric bifurcating solution branches) when dealing with 
the exact reduced map. 



8. Approximate breather solutions and comparison with numerical results 

In this section, the discrete breathers observed numerically in section [5] are shown to 
correspond to the different branches of solutions of the system of equations (I29ap - (l296ll 
discussed in the previous section. 

Assuming the case of a hard potential where /3 > in the potential ([2]), system 
(I29aj)-(!29l admits homocl inic solutions to with 6 = 0, and possesses homoclinic 
solutions to with a = Q ii k > 

64 

Homoclinic solutions of (]29ap - (l296l) with a = and b{^) = -\/2 Qq sech (5^) 
correspond to breather solutions with frequencies slightly above the optic band. The 
principal part of these solutions as the mass ratio m is close to | is provided by 
( 13T|) . Heavy and light masses oscillate at the same frequency 30;, and light mass 
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displacements are larger than heavy mass displacements (see figure [TJ bottom row). 
The existence of these solutions has been proved in [8] for an arbitrary mass ratio m in 
the small amplitude limit (for frequencies slightly above the optic band). Among the 
solution families obtained in [8], one family has the symmetry X_„ (t) = X„ (t) for even 
potentials, whereas the principal part of (13T1) only provides X_„ (t) — X„ (t) = O (e^). 
To better describe this symmetry, it is more convenient to rewrite (|3T|) in the form 



(t) = ~ v f '^0^ sech {6 en) cos (Stut) + O(e^), 

V (48) 

^2x1-1 (t) = aoe sech [6e [n - i)] cos (Scut) + O(e^), 

valid for fixed k > — ^ and e ~ (see figure [TJ bottom, for a comparison of this 
asymptotic expression with a numerically computed discrete breather solution). Fixing 
K and varying e is equivalent to moving parameters on the curve defined by 

- - = - - 2. (49) 
K \ 4/ m 

All these curves pass through (cj^, m) = (|, |). The first limiting curve F_i/64 is tangent 
at this point to the line Quj"^ = 2 (1 + m) (where 3u reaches the top of the optic band), 
and the second limiting curve Fqo corresponds to the line = 2m (top of the acoustic 
band) . 

Homoclinic solutions of fl29a|) - fl29l with 6 = and a (^ = ai sech (^^^j 
correspond to gap breather solutions with frequencies slightly above the acoustic 
band. Heavy and light masses oscillate with the same frequency u and heavy mass 
displacements are larger than light mass displacements (see figure [TJ top). The existence 
of such solutions has been proved in [8] in the small amplitude limit (for frequencies 
slightly above the acoustic band), and for mass ratio m outside intervals of resonance 
in which the acoustic band and an image of the optic band intersect (see shaded areas 
in figure [I]). As observed in section iGAl for an even interaction potential only the odd 
multiples of the optic band constitute forbidden zones. The mass ratio m = ^ precisely 
belongs to the boundary of one of these zones. System flHTl) formally extends the set of 
solutions obtained in [8] , since it provides approximate acoustic breather solutions with 
parameters close to {uj^,m) = (|, |), k > — ^ being fixed and e close to 0, 

(t) = (-1)" V2aiesech (^en) cos (cut) + 0{e'^), 
X2n-i (t) = 0(e2) 

(see figure [71 top, for a comparison of this asymptotic expression with a numerically 
computed discrete breather solution). 

Homoclinic solutions with two nonvanishing components a, b mix the two types of 
breather solutions (see figure [3, middle). According to fl3T|) . a period tripling bifurcation 
occurs when a solution with a nonvanishing a-component emerges from a solution 
consisting of a 6-component only (figure [71 from bottom to top). The case when a b- 
component emerges from a solution consisting of an a-component only does not produce 



(50) 
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a frequency change but it modifies the spatial pattern of hght mass oscillations and 
increases their amplitude. 

Let us now provide approximate expressions for the additional discrete breather 
solutions near the period tripling bifurcation. 

The bifurcating homoclinics defined by (I36l) - (l38ll and the approximation (!3Til 
provide approximate solutions 



(t) = —y^aoesech. {5en) cos (3ct;i(:) + (— l)"'aoeA sech* {6en) cos (co't) + 0(e^ + eA^), 

(51) 

^2n-i (t) = ^floesech [6e {n - i)] cos {3ut) + O (e^ + eA^) , 



parameterized by A ~ and e <^ A, with 

,2 



= — - 2, (52) 
m 

K-z^o = Q + 16fi:o)'A2 + 0(A^), (53) 



UJ 



'-\ + (54) 



4 



In figure [71 third row, we compare this asymptotic expression for e > and A > (right 
column) with a numerically computed discrete breather solution (middle column). The 
numerically computed discrete breather in the left column would correspond to e > 
and A < 0. Note that the map (A^,e^) h-^ (u;^,m) is locally invertible, and the value of 
(A, e) is determined by the choice of (cu^, m) specified in the figure. Moreover, bifurcating 
homoclinics with A 7^ locally exist for parameters above a curve m = mpT {00) passing 
through (u;^,m) = (i, |), with the local approximation (l40l) . 

The local bifurcation around acoustic breathers can be described in a similar way. 
The bifurcating homoclinics defined by (H^ - (H6|) and the approximation (|3T1) provide 
approximate solutions 

{t) = (— 1)" V^aie sech (^^en^ cos (cut) — -^OieAsech* (^^en^ cos {3ujt) 

+0 (e^ + eA^) , (55) 

X2n-i (t) = flieA sech" {n - |) cos (Siut) + O (e^ + eA^) , 

parameterized by A ^ and e <^ A, with 

= — - 2, (56) 



m 

K 



n, = J^^e\^ + 0[\'), (57) 
= i + Ke^ (58) 

(see figure [71 second row, for a comparison of this asymptotic expression with a 
numerically computed discrete breather solution, in the case e, A > 0). Moreover, 
bifurcating homoclinics with A 7^ locally exist for parameters below a curve m = 
'T^HM (i^) passing through {uj'^,m) = (^,|), with the local approximation (jlTj). At 
leading order in e and A, the frequency of light masses is three times the frequency u of 
heavy masses. This behavior corresponds to what is observed numerically for discrete 
breathers in this region of parameter space (see figure [TO] for an illustration). 
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Figure 10. Mass displacements as a function of time, plotted over one oscillation 
period for the five central particles of a gap breather close to an order three frequency 
resonance between the optic phonon with wave number k = and the acoustic phonon 
with wave number k = tt. The small masses (second and fourth curve from above) 
oscillate at leading order with three times the frequency of the heavy ones. Higher order 
corrections account for the slight modulation visible in the small masses' trajectories. 

Comparing the numerically computed discrete breathers with the analytic ap- 
proximations (H8l) . fl50l) . flSTl) . and fl55l) obtained by means of a center manifold 
reduction as plotted in figure [71 we found very good qualitative agreement of the 
results. Furthermore, excellent agreement is found when comparing the asymptotic 
behaviour of the bifurcation lines (l40l) and (l47j) of the reduced system (I29ap - (1296I) with 
the numerically computed bifurcation lines close to the "corner" formed by the phonon 
bands (see figure [9]). The quantitative agreement of the numerically computed discrete 
breathers with their analytic approximations depends, as expected, on how good the 
assumptions of the asymptotic analysis are met. For the weakly localized breathers 
in the upper row of figure U\ the oscillation amplitudes from numerical and analytic 
computations agree very well, whereas for the stronger localized breathers plotted in 
the bottom row the amplitudes differ by a factor of approximately |. Considering 
parameter values (cu^, m) closer to the phonon band edges, the localization strength of 
these breathers should decrease, leading to an improved agreement of numerical and 
analytic results. Unfortunately, we encountered problems with the convergence of the 
algorithm used for the numerical computation of the discrete breathers which prevented 
us from making such a comparison. 
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9. Summary and conclusions 

Discrete breathers in a diatomic Fermi-Pasta-Ulam chain have been investigated by 
means of numerical as well as analytical methods. The technique of numerical 
continuation, based on Newton's method for computing zeros of a function, was 
employed for the numerical analysis. Bifurcations between different types of discrete 
breathers were found, with the mass ratio m and the breather frequency u as bifurcation 
parameters. A common feature of the bifurcation lines observed is that they emerge 
from corners or edges formed by phonon bands and their images in the parameter plane 
{u"^, m). 

Above the optic band, a bifurcation between symmetric and asymmetric discrete 
breathers was investigated. The bifurcation line, starting from the upper edge of the 
optic band at mass ratio m = 1, was determined. 

In the gap between the acoustic and the optic band, discrete breathers are believed 
to exist only outside certain "images" of the optic band, i.e., when their frequencies 
are nonresonant with the optic phonons. These images cross the upper band edge of 
the acoustic band, forming "corners" in parameter space. We chose a particular corner 
at (ci;^,m) = (|, |) and investigated small amplitude discrete breathers in its vicinity. 
Two types of bifurcations, with bifurcation lines approaching the corner, were observed 
in this region. One is a period tripling bifurcation, the other is a local bifurcation 
around an acoustic breather increasing light masses oscillations with respect to heavy 
masses oscillations. Varying parameters across these two bifurcations, we have computed 
new types of breather solutions interpolating between the known optic and acoustic 
breathers. 

Using the center manifold reduction theorem, we have proved locally [for (cu^, m) ^ 
(i, I)] that all small amplitude (time-reversible) periodic oscillations are determined by 
a four-dimensional mapping. Homoclinic solutions of this reduced system have been 
studied formally by passing to a continuum limit, and a system of differential equations 
was derived, providing exact asymptotic expressions for discrete breathers close to 
(ci;^,m) = (|,|)- These expressions allow to identify the different types of breathers 
observed numerically with different branches of solutions of the set of equations. This 
analytic result confirms that, as suggested by the numerical data, the bifurcation lines 
emanate from (u;^, m) = (i, |), i. e., at the corner formed by the acoustic band and an 
image of the optic band. Moreover we have analytically computed the slope of the two 
bifurcation lines at this point. 

We want to conclude with remarks concerning the above results. 

Although numerical as well as analytic results were derived in the vicinity of a 
particular crossing of upper edge of the acoustic band with that of the second image of 
the optic band at (cj^, m) = (|, |) , similar bifurcations should take place close to all the 
similar "corners" in parameter space (with the harmonic cos (3t) in (!28|) replaced by an 
harmonic cos ((A; + 1) t) near the k-th image of the optic band). Moreover one should be 
able to detect 5 additional bifurcation curves near (cj^, m) = (|, |) , corresponding to the 
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additional critical values fii, . . . ,f25 leading to bifurcations around acoustic breathers 
(we have not explored this problem numerically). 

Lastly, an interesting result would be to prove for the exact map (l22^-(!22dll the 
existence of homoclinic orbits to close to the approximate homoclinics of section [71 
This would automatically imply by theorem [1] the existence of corresponding exact 
breather solutions for the FPU system. 
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Appendix 

In this appendix we report the details of the calculation of the approximate solutions of 
the four- dimensional recurrence relation fl22ap - fl22(ij) obtained by inserting the expansion 
( l25l) into these recurrence relations and identifying the different harmonics and powers 
of e. Bringing together different terms in (125|) we find 

a„ = i-ireA, (0 + e'A2 (0 + e'E, (0 , (1.1a) 

=e5.(0 + (-ire'{;e(0, (l-lb) 

cn = eC, (0 + (-1)" e'C^p + i-ir e'G, (0 , (l-lc) 

dr.+, = {-iTe'D, (0 + e'D^ (0 + e^H, (0 . {l.ld) 

The functions A, (0, (0, (0, (0, (0, F, (0, (0, (0 are of O (1) as 
e — s> 0. We have in addition 

B, (0 = -Sa (0 + O (6) . (1.2) 

In the following, we shall omit the e-dependency in the notation. Inserting Ansatz 
fll.lal) - p.ldl) into the recurrence relations fl22al) - (l22d[) . one obtains D2 = C2 = A2 = 0, 
= O (e^) and the equations 

= |/3e3 (-1)"+^ A (0 [A^ (0 + (-1)"^' A (0 {B (0 - C (0) + 2 (5 (0 - C {i)f] 

+ (-1)" 6^ {D (0 -D{i-e))+ (-1)" A (0 + 2e'H (0 + O(e^), (1.3a) 
= 2/?e3 [\ (-1)"+^ A^ (0 + (0 {B ii) - C (0) + I (5 (0 - C {i)f] 

+e {B iO + B (e-6)) + (96^ + 16) eC (0 + 16 (-1)" e'G (0+ O(e^), (1.3&) 
= I (-1)"+^ e'D (0 - le'H (0 + (-1)"+^ e (A (^ + e) - A (0) 

+2e'E (0 - iPe'A' (0 (B (0 - C (0) + O(e^), (1.3c) 
= {9Ke' + i) eB (0 + e (C (^ + e) + C (0) 

-3Pe' [A' (0 {B (0 - C (0) + i (5 (0 - C iOf] + O(e^). {1.3d) 
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Using (\1.3b\i we get 

G (0 = ^/5^' (0 + 0(6). (1.4) 

Now fll.36p allows to eliminate C as a function of A, B up to order e^. Substituting the 
result into (ll.3(ip yields 

= 5(^ + 6) -25 (0 + i?(^-e)-9e2(i + 16«:)i?(0 

+ 243^^2 ^2^2 ^ ^ 81^3 (^)] + 0(6^). (1.5) 

Similarly, fll.3ap and fll.Scp lead to 

H{0 =-!/?^^(0(i?(0-C(0) + O(e) 

= -1/3^2(05(0 + 0(6), (1.6) 
E (0 = (0 + \PA' (0 (5 (0 - C (0) + 0(e) 

= ^/3A2 (05(0 + 0(e), (1.7) 

e5(0 =-|(A(e + e)-A(0) + O(e^). (1.8) 
Using ([L8]) and ([Ol) in (Ool) yields 

= A (e + e)-2A (0 + A (e - e)-|e2A (0 + f /5e' [A^ (0 + (0 (0] +0(e=^).(1.9) 

We proceed by expanding A, B, C, D, E, G, H up to first order in e, while ^ is kept 
fix. Recalling that A (0 = Aq (0 + O (e), B (0 = Bq (0 + O (e) and using the same 
notation for the other variables, one obtains from the above equations as e ^ 

= ^ - 1^0 + f /5 {Al + %AoB^,) , (1.10a) 

= ^ - 9 (i + 16k) B, + {2AlB, + f^B',) , (1.106) 

Co = -|5o, (1.10c) 

Do = -1^, (l.lOrf) 
d^ 

Eo = ^PAlBo, (l.lOe) 

Go = ^,PAl, (1.10/) 

Ho = -%PAlBo. il.lOg) 
Going back to fll.lap - Ql.ldp . we end up with 

an= (-l)"eAo(en) + 0(e2), (1.11a) 

6„ =e5o(en) + 0(e2), (1.116) 

On = -|e5o(er2) + 0(e2), (1.11c) 

cin = f(-l)"62^(en) + 0(e3), (l.llrf) 

where Ao,Bo are solutions of system fll.lOap - (]1.106p . In this way we have obtained a 
formal continuum limit of the four-dimensional recurrence relation fl22ap - (122dp which 
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provides approximate solutions valid for 77 0. 
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